Zener double exchange from local valence fluctuations in magnetite 



o 
o 



u 

a 

O 

o 



> 
in 

<N 
(N 

o 
o 



R. J. McQueeney/'^ M. Yethiraj,3'E S. Chang,^ W. Montfrooij T. G. Perring,^ J. M. Honig,^ and P. Metcalf^ 

^Department of Physics & Astronomy, Iowa State University, Ames, Iowa 50011 USA 
^Ames Laboratory, Ames, Iowa 50011 USA 
^Center for Neutron Scattering, Oak Ridge National Laboratory, Oak Ridge, TN 37831 USA 
^Department of Physics and Missouri Research Reactor, 



University of Missouri, Columbia, MO 65211 USA 
^ ISIS Facility, Rutherford- Applet on Laboratory, Clarendon, 



UK 



^Department of Chemistry, Purdue University, West Lafayette, IN Ji.1901 USA 

(Dated: July 15, 2007) 

Magnetite (Fe3 04) is a mixed valent system where electronic conductivity occurs on the B-site 
(octahedral) iron sublattice of the spinel structure. Below Ty — 122 K, a metal-insulator transition 
occurs which is argued to arise from the charge ordering of 2+ and 3+ iron valences on the B- 
sites (Verwey transition). Inelastic neutron scattering measurements show that optical spin waves 
propagating on the B-site sublattice ('-^80 meV) are shifted upwards in energy above Ty due to the 
occurrence of B-B ferromagnetic double exchange in the mixed valent metallic phase. The double 
exchange interaction affects only spin waves of A5 symmetry, not all modes, indicating that valence 
fluctuations are slow and the double exchange is constrained by electron correlations above Ty . 

PACS numbers: 71.30.+h, 75.30.Ds, 75.30.Et, 78.70.Nx 



Magnetite (Fe304) is the prototypical example of a 
metal-insulator transition with a charge-ordered (CO) 
insulating ground state (Verwey transition). Since its 
discovery nearly 70 years ago,[l] the driving forces be- 
hind the Verwey transition are still not completely 
understood. [3, 0] Magnetite has a cubic inverse spinel 
crystal structure containing two different symmetry iron 
sites; the A-site resides in tetrahedrally coordinated oxy- 
gen interstices and has stable valence (3d^, Fe^+), the two 
B-sites have octahedral coordination and a fractional av- 
erage valence of 2.5+. The ferrimagnetic structure con- 
sists of ferromagnetic A- and B-sublattices aligned an- 
tiparallel to each other (Tc = 858 K). Below Ty = 122 
K, magnetite undergoes a metal-insulator transition re- 
sulting in a decrease of the conductivity by two orders- 
of- magnitude. The model that has persisted over time 
is that extra electrons forming Fe^+ ions (3d^) hop to 
neighboring Fe^+ sites on the tetrahedral B-sublattice 
network and give rise to electrical conductivity. Ander- 
son argued that short-range ordering of Fe^+ and Fe^+ 
exists above Ty due to significant intersite Coulomb re- 
pulsion and frustration on the B-site sublattice. [J] The 
short-ranged electron correlations maintain local charge 
"neutrality" (2.5+ average valence on each tetrahedron), 
thereby restricting charge hopping and conductivity. [5| 
Eventually, Coulomb repulsions win out at low tempera- 
tures, resulting in long-range CO of Fe^+ and Fe^ [6| in 
a process reminiscent of Wigner crystallization. [7] Weak 
elastic and orbital interactions [i| induce monoclinic lat- 
tice distortions^whose complexity obscures the details of 
the CO state. [9] Even the validity of the CO model has 
been questioned. [10] However, recent neutron fn\ and 
resonant x-ray scattering measurements [12] appear to 
converge on fractional CO. 

In this Letter, we provide strong evidence for Ander- 



son's original picture of the Verwey transition in mag- 
netite as arising out of short-range electron correlations 
in the mixed valent (MV) phase. Valence fluctuations 
occurring on the B-sublattice modify the magnetic ex- 
change and affect spin waves propagating on Fe B-sites. 
Inelastic neutron scattering measurements of spin waves 
reveal that B-site optical spin waves are shifted up in 
energy and broadened above Ty due to ferromagnetic 
double exchange (DE). Ferromagnetic DE arises from 
real charge transfer processes in MV materials, a good 
example being the ferromagnetic metallic state in the 
manganites.[13] For fast electron hopping in the band 
limit, the average DE should uniformly affect all B-B 
pairs. However, our results show that only spin waves of 
a particular symmetry are affected by DE, implying the 
presence of Zener DE [3] arising from slow (< 80 meV) 
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FIG. 1: Spin wave dispersion of magnetite above Ty. Black 
(gray) symbols are inelastic neutron scattering data from 
MAPS (HB-3). Lines are results from a Heisenberg model 
with parameters discussed in the main text. The hatched area 
contains very broad B-site spin waves. 
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FIG. 2: (a) Neutron intensity in the [HKO] plane at 130 K. The color scale indicates neutron intensity (red the highest). Thin 
lines show the Brillouin zone boundaries of the cubic phase. Dashed lines show 0, 80, and 120 meV energy transfers (from left 
to right). The purple circled region indicates extremely broad and weak A5 symmetry optical spin waves at 80 meV. (b) The 
110 K data set below Ty. In the purple circled region, the A5 mode is more pronounced, (c) Calculation of the MAPS spin 
wave spectra using a Heisenberg model with 7^^=0.44 meV. Data and calculations have been smoothed by a boxcar averaging 
procedure, (d) Symmetry of B-site optical spin waves in different Brillouin zones at 80 meV; blue (A5), green (A2(B)), purple 
(V). Energy cuts shown in fig. 3 are indicated by arrows. 



electron hopping and restricted by short-ranged charge 
correlations. The local valence fluctuation symmetry can 
be inferred from the data and mirrors the eventual long- 
range symmetry of the CO state. [15] Inelastic neutron 
scattering measurements were performed on the MAPS 
instrument at the ISIS facility at Rutherford-Appleton 
Laboratory on a single-crystal of Fe304 weighing ~10 
grams. Details of sample preparation and characteriza- 
tion are given elsewhere. [16] The sample was mounted 
with cubic [HKO] as the primary scattering plane and 
data were collected with an incident neutron energy of 
160 meV and the incident beam at an angle of 25^ from 
the [110] axis. Time-of-flight neutron spectra were col- 
lected at 110 K (T < Tv) and 130 K (T > Ty) and scat- 
tered intensities were histogrammed into energy transfer 
{hco) and momentum transfer (HQ) bins. Data were sub- 
sequently analyzed using MSLICE and TOBYFIT 
[18] computer programs. The measured neutron scatter- 
ing spectra from MAPS (and also the HB-3 spectrom- 
eter at the High Flux Isotope Reactor at Oak Ridge 
National Laboratory) were used to determine the spin 
wave dispersion in magnetite. In the cubic spinel phase 
above Ty, the primitive rhombohedral unit cell con- 
tains six iron atoms (2 A and 4B), leading to six spin 
wave branches. Figure [T] shows the spin wave disper- 
sion along [100] as determined from Gaussian fits to spin 
wave modes. The dispersion is well represented by a 
Heisenberg model where antiferromagnetic AB superex- 
change, via oxygen, is dominant {Jab = — 4.8 meV) and 
is responsible for ferrimagnetism. Ferromagnetic BB ex- 
change arises from a combination of superexchange, DE, 
and direct exchange and is an order-of-magnitude weaker 
{Jbb = 0.69 meV). Weak antiferromagnetic nearest- 



and next-nearest-neighbor AA exchange is also present 
i'^AA = -0-35 meV, jjfj = -0.2 meV). The dispersion 
calculated from a Heisenberg model with exchange pa- 
rameters above and spins Sa = 2.5 and Sb = 2.25 is 
also shown in fig. [TJ Branch symmetries were identi- 
fied from the model spin wave eigenvectors and are given 
the following labels and descriptions; Ai (acoustic and 
steeply dispersing optic), A2(A) (optic spin wave on A- 
sublattice), ^2(3) (optic spin wave on B-sublattice), and 
A5 (doubly-degenerate optic spin waves on B-sublattice). 
Modes of A5 and ^2(8) symmetry at ~ 80 meV (in 
the hatched region of fig. [1]) propagate solely on the B- 
sublattice and were observed to be very broad and weak, 
as discussed below. Figures 2(a) and (b) show images of 
scattering intensity in the [HKO] plane at 130 K and 110 
K, respectively. The arcs and rings in the images cor- 
respond to the intersection of the neutron measurement 
surface and the dispersion surfaces. Figure 2(c) shows 
a calculation of the MAPS data using the Heisenberg 
model described above. Calculations simulate the sam- 
pling of spin waves in (Q,a;) space performed by MAPS. 
Features of the calculated and measured spectra are in 
good agreement (as anticipated from fig. 1) with one no- 
table exception. As indicated in Fig. 2(a), the measured 
B-site optical spin wave branch at ~80 meV is poorly 
defined in the region between [220] and [420] at 130 K. 
Comparison with the 110 K spectrum illustrates two fea- 
tures; (1) the B-site modes sharpen up below Ty in bet- 
ter agreement with the cubic model (despite the lower- 
ing of crystallographic symmetry), (2) other spin wave 
branches are essentially unaffected, indicating that Jab 
and Jaa are insensitive to the Verwey transition. These 
observations are supported by energy cuts through the 
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FIG. 3: Inelastic neutron scattering intensity as a function of 
energy transfer for cuts along; (a) [i^20], (b) [i^30], and (c) 
[H40] . Blue (red) symbols indicate measurements done at T = 
110 K (130 K). Blue (red) lines are Gaussian fits to the 110 K 
(130 K) data. The dashed Une is the background estimated 
from cuts along [if 60]. Panels (d) - (f) show the identical 
cuts as (a) - (c) as calculated from a Heisenberg model for 
magnetite with Jbb — 0.44 meV (blue) and Jbb — 0.69 meV 
(red). Spin wave modes are labelled by symmetry. 



MAPS data at different values of K and different tem- 
peratures, as shown in fig. 3. To improve statistics, cuts 
are obtained by averaging neutron intensity over a range 
of reciprocal space. Fig. 3 shows the following cuts made 
in the range -0.25 <L< 0.25; [if20] (1.75 <K < 2.25), 
[if30] (2.75 <K < 3.25), and [if40] (3.75 <K < 4.25). 
These reciprocal space cuts allow for specific symmetry 
selection of the B-site optical spin waves, as shown in fig. 
2(d). In the constrained (Q, ficj) measurement surface 
of MAPS, if is a function of hw^ as can be ascertained 
from the comparison of figs. 2 and 3. The energy resolu- 
tion full-width-at-half-maximum (FWHM) is calculated 
to be 3-4 meV for dispersionless modes in the energy 
range shown (using the TOBYFIT program). In all cuts, 
spin wave modes at 60 and 117 meV are independent 
of temperature and resolution limited, while ~80 meV 
modes are broader than the resolution. The energy cut 
along the [if 20] direction in fig. 3(a) shows that the A5 
symmetry branch at ^80 meV shifts up in energy and 
becomes severely broadened above Ty, as also apparent 
in fig 2. Gaussian fits to the data indicate that the mode 
energy shifts from 79.6(5) meV to 84(1) meV and the 
FWHM broadens by 60% from 17 to 28 meV above Ty. 
The B-site spin waves along [if 30] in fig. 3(b) (A'2(5) 
symmetry) and [if 40] in fig. 3(c) {V symmetry [19]) also 
have broad lineshapes with FWHM of 14 meV and 10 
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FIG. 4: (a) Neutron structure factor for B-site optical spin 
waves summed from 75-90 meV. Red (blue) circles are data 
at 130 K (110 K). Red (blue) lines are calculated from a 
Heisenberg model with Jbb =0.69 (0.44) meV. Horizontal 
bars indicate dominant symmetry, (b) Structure factor for 
A-site optical spin waves from 115 -120 meV. Black line is 
the result of Heisenberg model. 



meV, respectively, but relatively little temperature de- 
pendence. We also note the weak peak near 70 meV in 
figs. 3(b) and 3(c) is likely to be a phonon, based on 
examination of higher-Q data. The anomalous B-site 
spin waves can be explained by consideration of ferro- 
magnetic DE. The conduction electron associated with 
the Fe^+ ion is forced to remain oppositely aligned to 
the Fe^+ 3d^ core spin by the Pauli exclusion princi- 
ple (effectively Jnund — ^ 00). Since the B-sublattice 
is already ferromagnetically aligned due to strong Jab, 
spin polarized conduction results in DE[23| and leads to 
an increase in Jbb- In turn, the energy of B-site opti- 
cal spin waves increases. The enhancement of DE above 
Ty can be estimated from the energy of the dispersion- 
less A5 mode, 6JabSb + SJbbSb- Since Jab does not 
change, [16] the shift of 4.4 meV implies that Jbb in- 
creases from 0.44 meV to 0.69 meV. The energy cuts are 
compared to the Heisenberg model structure factor cal- 
culation in figs. 3(d)-(f) with Jbb = 0.44 meV and 0.69 
meV. Comparison to calculation shows that; (1) all B- 
site spin wave peaks are much weaker and broader than 
Heisenberg model predictions, (2) changing Jbb in the 
model shifts the energy of all mode symmetries, in dis- 
agreement with the observation that only the A5 mode 
has a large temperature dependence. Structure factors in 
the energy range 75 < hco < 90 meV and 115 < hco < 120 
meV and L-range —0.25 < L < 0.25 are shown in fig. 4 as 
a function of K. Fig. 4(a) further illustrates the anoma- 
lous A5 mode, which shows a minimum in the measured 
structure factor whereas the calculation predicts a max- 
imum. Fig. 4(b) indicates that A-site optical spin waves 
have little temperature dependence and agree quite well 
with Heisenberg calculations. It is clear that stiffening 
and extreme broadening of A5 modes above Ty cannot be 
explained by uniformly enhanced BB coupling as might 
be expected from a band approach. The failure of a band 
model is reiterated by estimating the DE arising from 



metallic conductivity, Jde ^ ^ t/2V2Sl [m, m where 
n ^ 1/2 is the concentration of charge carriers in the con- 
duction band and t the electron hopping integral. The 
effective one-band hopping integral from LSDA calcula- 
tion is ~ 130 meV[23], leading to an estimate of Jde ^ ^ 
meV that is significantly larger than the observed value 
{Jde - Jbb{T > Tv) - Jbb(T < Tv) = 0.25 meV). 
This disagreement emphasizes the importance of charge 
correlation and/or polaronic effects above Ty that lead 
to activated conduction and a reduction of the effective 
hopping integral. This picture is supported by optical 
conductivity measurements showing that the ~100 meV 
electronic energy gap below Ty (which suppresses DE) 
is replaced by a pseudogap above Ty.[2J, [25,] In this 
limit, the Zener DE picture of slow and local charge 
hopping between distinct valence states above Ty is ap- 
propriate and the closing of the electronic gap allows 
the coupling of spin waves to valence fluctuations lead- 
ing to broadened spinwave lineshapes. The predominant 
A5 spin wave symmetry would then imply that valence 
fluctuations have the same symmetry and occur along 
the quasi-one-dimensional [110] B-site chains, as illus- 
trated in fig. 5. The spin wave changes the alignment 
of neighboring spins, thereby impeding valence fluctua- 
tions which favor parallel spins. The presence of slow A5 
valence fluctuations coupling to the lattice has also been 
observed with neutron diffuse scattering and occur well 
above Ty [2^. In the MV phase, these local symmetry- 
breaking charge correlations are short-ranged. While we 
interpret our data based on the symmetry assignments 
for modes along the [100] direction, broadening every- 
where in the Brillouin zone indicates that the coupled 
B-site charge and spin fluctuations are local. The result- 
ing picture of valence fluctuations constrained by short- 
range ordering of Fe valences is similar to that proposed 
originally by Anderson. [4] In addition, the implication of 
slow electron hopping (<80 meV) supports more detailed 
theories of magnetite that invoke polaronic behavior for 
the Fe^+ charge carriers. Q Spin-charge coupling is large 
and may contribute to the polaronic binding energy. Fi- 
nally, the A5 symmetry of charge correlations are similar 
to the eventual frozen pattern of the CO state, which 
has been argued to be a combination of A5 X3 cubic 
representations. [15] 
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FIG. 5: Schematic drawing of (a) spin and (b) charge modes 
on the B-sublattice with A5 symmetry. In (b), the heavy 
arrow indicates direction of charge fluctuations. 
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